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NATIONAL ADVISORY COMMITTEE FOE AERONAUTICS 


TECHNICAL NOTE NO. i486 


STRESS DISTRIBUTION IN A BEAM OF OSTEOTROPIC MATERIAL 
SUBJECTED TO A CONCENTRATED LOAD 
By c. B. Smith and A. V. Voss 


SUMMARY 


Mathematical expressions have teen derived for the stress distribu- 
tion in a wood beam of rectangular cross section subjected to a concentrated 
load. The orthotropic nature of wood was taken into account in the deriva- 
tion. The stress distributions were expressed in terms of infinite Beries. 

A method of reducing the infinite series to the sum of a finite series and 
a closed form was described. 

The mathematical ly determined distribution of horizontal shear in 
the vicinity of a concentrated load was compared with the actual distribu- 
tion obtained by a test of a Sitka spruce beam of rectangular croBB section. 


INTRODUCTION 


A number of experiments on the bending of wood beams conducted at the 
Forest Products Laboratory have shown some results that are not explainable 
by the elementary theory of bending. It is well known that the state of 
stress that is produced in the interior of a beam* slightly bent by any 
forces, may be approximated by the elementary theory of bending at all points 
that are at a considerably large distance from any place of loading or of 
support. But the stress distribution near a concentrated load or a place 
of support is not easily determined. Hence, in applying the usual beam 
theory to wood beams, discrepancies occur ’in the neighborhood of concen- 
trated loads. This analy sis is an attempt to present for wood beams a more 
nearly e xa ct mathematical derivation of the stress distribution near a 
concentrated load that is obtained from the elementary theory of bending, 
in order to explain some of the discrepancies that may arise in the bending 
of wood beams . 

This work was conducted at the Forest Products Laboratory under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 
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MATHEMATICAL ANALYSIS 
Method of Analysis 


The beam discussed is assumed to be an orthotropic solid in the form 
of a long, thin, rectangular plate having its edges parallel to two perpen- 
dicular axes of elastic symmetry lying in the plane of the plate. In the 
analysis wood is considered to be orthotropic. (See references 1, 2, and 3«) 
For mathematical simplicity the thickness of the beam is assumed small as 
compared with the vertical depth of the beam so that the problem can be 
treated as one of plane stress. 

The beam is taken ter be infinitely long and. to be subjected to a 
periodic normal load on the upper and lower faces (references 4 and 5) • 

The results obtained are then extended to various types of loading and end 
conditions for a beam of finite length. The formulas obtained apply to 
beams of any thickness that is small in comparison with the depth If the 
load Is considered to be given per unit thickness. 

The x-axis Is taken along the middle line of the beam, and the equation 
of the upper and lower faces of the beam Is taken to be y = ± h, as shown 
by figure 1. 


Sinusoidal Loading of an Infinite Beam 


For the state of plane stress in the orthotropic beam, the stress 
function is given by a suitable solution of the following differential 
equation. (See reference 6.) 

fih; + 2 k aV - + =, o (i) 

«x* 8x 2 8„ 2 ^ 


where 


and 



T} = ey 


( 2 ) 

(3) 


6 



(4) 


In equation (.3) 
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Also, 

Ey moduli or elasticity in x— and y-directions, respectively 

jjLjy. modulus of rigidity associated with xy— plane 

CT^y Poisson 1 s ratio associated with, stress in x-direction and 

strains in x— and y-directions 

For the components of stress the notation Yy, and Xy 1 b used 

as in reference 7* First the team is subjected to the surface force 




cos mx 


*y - 


on the edges y = 1 h. 


A solution of equation (l) that can he made to satisfy these houndary con- 
ditions is 


F = ^Aj cosh hkitj + Bj_ cosh mflT^ cos mi 


where 


a = W k + \[k 2 — 1 


P =\|k K 2 — 1 


and it is important to note that a£ = 1. 
The resulting stress components are 


= Yy = -m 2 ^A-j_ cosh mat) + Bj cosh m0T^ 


cos mx 


€ 2 s X x =■ e^m 2 ( A^a. 2 cosh mail + BjP 2 cosh mBr[\ cos mx (9; 

6t| 2 V ' 

— c - = Xy = cm 2 ^A^a sinh mat} + BxP sinh mBt)^ sin mx. (10) 


From equation (5)* the value of 2^. on the surface leads to 
A^a sinh mcc-eh + B^P sinh mfkh = 0 


B ! = 


Aqa sinh maeh 
p sinh mBsh 
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Also from equation ( 5 ) j the value of Yy on the surface leads ter 

^A-^ cosh maeh + B^ cosh mPeh^ cos mx = gBT'cos mx (12) 

On solving equations (ll) and ( 12 ) for A^_ and B-j_, it is found that— 

HP sinh mJ 3 eh 


An = 


■1 T~. 

2m . (a sinh maeh cosh mpeh — p sinh mpeh cosh maeh) 


( 13 ) 


and 


B, = 


-Ha sinh maeh 


^ 2m2 (a sinh maeh cosh mPeh — p sinh mpeh cosh maeh) (l^-) 

The state of stress is now given by 

Y y = — ^p sinh mPeh cosh marj — a sinh maeh cosh mPr^ cos mx (15 ) 
^a, sinh mPeh cosh marj — p sinh maeh cosh mP^ cob mx (l6) 
He 


x x = *£ 
X 2Dl 


Xy 2 Di 


^sinh mPeh sinh mat) — sinh maeh sinh mPi)^ sin mx (17) 


where 


= a sinh maeh cosh mPeh — P sinh mPeh cosh maeh 
Next the boundary conditions are taken to be 

Yy = cos mx (j = h^ 


Y = — =kE" cos mx 

y 2 

Xy = 0 


(y = - h) 
(y = ± h) 


(18) 


A solution of equation (l) that-can be made to satisfy these condi- 
tions is • 


F = ^2 sinh mar) + B 2 sinh iuPt]^ cos mx 
and the resulting stress components are 


6% 


8x 


r 2 / N 

~2 = Yy = — m f A 2 sinh marj + B 2 sinh mPt] j cos 


mx 


( 19 ) 


( 20 ) 
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g 2 _ e 2 m 2 ^^0,2 git^h nua,^ + b 2 P 2 sinh mPr)^ cos ms 

— 6 5-JL = Xy = em^ ^A^a. cosh mar) + B 2 P cosh mPr|^ sin mx 


cos no: 


From, equations (l 8 ) and (22) it results that 

AgOt, cosh maeh + BgP cosh mpeh = 0 


b 2 = - 


Aoa cosh maeh 


P cosh mpeh 


Also from equations (l 8 ) and (20), it follows that 


n? ^A 2 sinh maeh + B 2 sinh mPeh^ = — iff 


On solving equations ( 23 ) and (24) for Ag and B 2j it is found that 


■^2 " 


— HP cosh mP 6 h 

2 m 2 (p sinh maeh cosh mPeh — a cosh maeh sinh mpeh) 


_ Ha cosh maeh r^\ 

b 2 = — ( 26 ) 

2 m (P sinh maeh cosh mPeh — a cosh ma 6 h sinh mpeh) 

By substituting from equations ( 25 ) and ( 26 ), the stress components 
become 

Y = f p cosh mPeh sinh mai) — a cosh maeh sinh hiPtA cos mx (27) 

y 2D 2 J 

2_ = — f a cosh mpeh sinh man — P cosh maeh s'inh mPT^ cos mx 

“2 V ) (2g) 

2^ = — ^cosh mpeh cosh mar) — cosh maeh cosh mPr)^ sin mx ( 29 ) 


where 

D 2 = P sinh maeh cosh mpeh — a sinh mpeh cosh maeh 

For isotropic material, it follows from equations (2), (4), and ( 7 ) 
that a, P, and e are equal to unity. For this value of a, p, and e, the 
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stresses given in equations ( 15 ) to ( 17 ) and ( 27 ) to ( 29 ) reduce to the 
ones that have teen found for the isotropic case. (See references 4 and 5*) 

The two problems just solved give periodic sinusoidal load distribu- 
tions along the upper and lower faces of the beam. It is now evident that 
by means of Fourier series the beam can be subjected to a very general 
type of loading. 


Concentrated Loading of - an Infinite Beam 


Consider first the infinitely long beam acted on by equal and similarly 
directed loads distributed over equal intervals of length a. This distri- 
bution of loads may be regarded as a continuous load of the type 

^r(x) = ... ft (x — 2a) + 0(x — a) + ft>(x) + 0(x + a) + 0(x + 2a) + ... ( 30 ) 

where fi(x) id a suitably restricted even function of x. It follows that ~ 
t]t(x) is an even function of period 2a and that ty(x + a) = tKx — a) = ijf(x) ■ 
These requirements are satisfied by the type of load used in reference 5, 
in which the definition of jb(x) was chosen as 


0(z) 


5 1 


( 31 ) 


This expression represents a unit load, since the area between the curve 
and the x-axis is unity. For small values of & It is suitable for the 
approximate representation of a load applied over a small area by a curved 
loading block. The components of stress associated with a point load are 
obtained as the limit, as 5 approaches zero, of the expressions for these 
components in terms of the parameter &. These limiting expressions repre- 
sent the exact solution of the problem of determining the stress distri- 
bution associated with a point load. For finite but small values of & 
the expressions for the components of stress give the approximate distri- 
bution of stress associated with a load distributed over a small area. 


The expression for ijf(x) (equation ( 30 ) ) can be represented by a 
Fourier series of the form 


where 


t(x) = bQ + \>2 cos + bij. cos + ... 


a ° ■ c 

2 r 


\jr(x) dx 

K = +(*> cos ES 


t 'l£> 


dx 


(r = 2, 4, 6 . . 


> 


(32) 


( 33 ) 
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Substituting equation ( 30 ) la the first of these integrals gives 

f pa pa pa 

b Q = ^ . . . + J 0(x — 2a) dx + J 0(x — a) dx + J 0(x) dx 


pa na *| 

0(x + a) dx + J 0(x + 2a) dx + . . .1 


(3b) 


Let 


x — 2a = x-Q 
x - a = x_ 2 _ 

x = Xq 

x + a = x n 


> 


(35) 


Then 


^0 = i 


[ p-a AO na 

... +J ^ 0 (x_e) dx_g + J 0 (x_ 1 ) dx_-L + / 0 (xo) dxo 

02a H3a -j 

+ J 0 (xx) dx x + J 0 (x s ) dxg + . . . 


or 


*0 - 41 


(36) 


By similarly substituting equation (30) in the second of the integrals 
of equation ( 33 ) } it follows that 


[•' *1 

t 


0(x — 2a) cos dx +J* 0(x — a) cos • dx 


ia pa 

+ / 0(x) cos ESS. dx + J 0(x + a) cos ^SS. dx 


+ / 0(x + 2a) cos dx + 


-] 


where r = 2, h, 6 .... 
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Again make the change of variable indicated in equation (35)* It 
results that 




+ f 0(x_ 2 ) COS ^ (x _ 2 + 2a) 
U — 2 a 


J 0 (x__ 1 ) cos Ef (x _ 2 + a) dx _ 1 +f 0 (xq) 008 ^ *0 ^ 


fl2a 


4 / 0{x;q) cos ~ (x^ - a) cLxt|_ 4 / 0(x 2 ) cos ~ (x 2 — 2 a) dx£ 


U0 

P3a 


l/a 

where r = 2 , 4, 6 . . 
or 


U 2 a 




h-.. = - 


0 (x) cos — dx 


(37) 


where r = 2 , 4, 6 .... 

On substituting from equation (31), equations ( 36 ) and ( 37 ) become 


*0 " a 


rrt 8 


*r - f e 


2 . — ar 


(33) 


where r= 2 , 4, 6 .... 


Next consider the infinitely long beam acted upon only by a continuous 
load of the type 

i$-(x) = ... + 0(x — 2a) — 0(x — a) 4 - 0(x) — 0(x 4 a) 4 0(x 4 2a) ... (39) 

where as before 0(x) is a suitably restricted- even function of x. The 
function Wx) is an even function of period 2a. Further, 
ty(x 4 a) = ij r (x — a) = -Hjr(x). Consequently, ^(x) can be represented by 
the Fourier series 

tJt(x) = b]_ cos ~ 4 b^ cos + b^ cos ^ 4 ... (40) 

where 

t s = | t(x) cos —• dx (s-l, 3 , 5 . .A (4l) 
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On substituting from equation (40) , equation (4l) becomes 

0(x — 2a) cos S|£. dx — -P 0(x — a) cos dx 


*■ ■ ■ [" •£ 


vf 


-r 


+ j 0(x) cos — — dx — J 0(x + a) cos —jp dx h f 0(x + 2a) cos dr 


{ 


where s = 1, 3, 5 • • • • 

By making use of equation ( 3 % it follows that 


a a 


fl— a 


+ J 0(x-e) cog ” (x_g + 2a) dx_ 2 


sitxo . 

COB -g— dx Q 


-J cos ^ (x^ + a) dx_^ +£ 0(xq) 

p 2a A 3a ~l 

—J 0 (^2^) cos ^ (x-j^ — a)dr^ +J 0(xg) cos (x.^ — 2 a)dX 2 •* J 


where s = 1, 3 3 5 ■ . . 
or 

i 

where s = 1, 3, 5 • • • • 


-iS. 


0(x) cos ^S. dx 


00 


By making use of expression ( 31 ) for 0(x) , it follows that 

sit5 

b — — e~ a 
D s - a ® 


(42) 


where s = 1, 3 } 5 • • • and 6 is taken very small in. the applications that 
follow. 


Beam of Finite Length 

It is now possible to investigate the main problems of this paper, 
that is, finite beams loaded in various ways and having either clamped or 
freely supported ends. 

Take the load to consist of a series of isolated loads -W on the 
upper face of the beam, with a series of equal and opposite supporting 
pressures on the lower face halfway between the loads, as shown in 
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figure 2 (a) . Tills distribution of forces may be resolved Into the sum of 
the two distributions shown in figures 2(b) end 2(c) , which are most con- 
veniently treated separately. 

In order to calculate the effect of the forces shown in figure 2 (c), 

it follows on reference to equations (15) j (30)* and (38) that m = — , * 

s. 

and it is necessary to write 


or 


H 


ritS 

H = — — e ~ 
a 


( 43 ) 


accordingly as r = 0 or r > 0 and sum for even values of- r. Thus, 
the stress components are, in figure 2(c); 


T -2 
y a 


-i + 


ritS ^ 

£31 h silfc 00( * = 

Hl V a 


_ 3 C =2 J • • • 

-a sinh EJfosfe cosh cos I 

a a / a I 


.. 2 , “ ritb /■ 

_ K£_ y e~~ (a sinh &&£& COB h ESSU 

a - 4 — m \ a a 

r=2, ... 1 


— 3 sinh 


coah £2 &lS CO b £5* 

a a J a 

rrtS 


J. y> 


r=2, . . . 


— sinh sinh 


rjrgjQ.^ 


sin 


rrcs 




( 44 ) 


where r has even values and 

Dj_ = a sinh 2^2^ cosh EUSik _ p Blnh £*21^ oosh Z 2 *sk 


( 45 ) 


In order to calculate the effect of the system of forces indicated in 
figure 2(b), it appears on reference to equations (27), ( 40 ), and ( 42 ) that 
it is necessary to write 

srtS 


a 



NACA TN Wo. i486 


11 


and sum for odd values of s. Thus, the stress components are, for the 
case shown In figure 2(h) : 


Y = -* 
7 a 


S — 1 , • • a 


srtS f 

e a (p cosh sinh S X& H 

Do '• a a 


_ a C osh sinh ^ 0 os && 

a a / a 


0 


W6 2 


Srt6 


s=l, ... 


(■ 


e a {a cos h sinh 

Dr ' a 


_ p cosll Sitaeh slnh saga 
a a 


COB 


*y“ T- 


sjtS 

a 


s=l, ■ 


/ , sitpeh . sjraTi 

(cosh — - — cosh L 

a a 


- cosh C osh Sf^sin SS 


a 


a 




where s has odd values and 

Do = p amt S22£k cosh _ a slnll cosh ^ 


a 


I 


(46) 


J 


( 4 ?) 


If now S is taken to he very small, the sum of the two stress dis- 
tributions (equations (44) and (46)) would give a stress distribution that 
very closely approximates one arising from a series of loads ¥ distrib- 
uted over small areas, as shewn in figure 2(a). It Is evident from 
figure 2(c) that the stress components given by equation (44) will be 
relatively unimportant except in the immediate neighborhood of the place 
of application of the forces. 

It is now possible to draw some conclusions concerning the flexure of 
finite beams under concentrated loads. Consider the portion represented 
by OS in figure 3» It closely approximates a beam of length 2a clamped 
horizontally at the ends and carrying a load ¥ at the center, 

* 

Again, consider the portion EE. This part closely approximates a 
beam of length a supported by vertical shearing forces of amount 

on the two terminal sections, having zero bending moments at P end. R 



H 

and. carrying a load W at the center. Thus section ER represents approximately a beam 10 

of length a, simply supported at the ends and carrying a load V at the center. 


A Simple Method for Computing Equations (44) and (46) 


It is possible to express approximately each stress component given in equations (44) 
and (46) in two parts, each having a finite number of terms. The approximation can be made 
as close as desired. As an illustration of the method of transforming the expressions for 
the stress components in this manner, the process is carried through in detail for the stress 

component ly given in equation ( 44) . Ey dividing both numerator and denominator by 

cosh 32S4 — cosh c \ it results that the stress component X_ in equation (44) becomes 

Q. 6L J 



nt6 

a 


f ^-\slnh ^Sai vsinh 

(tanh I— - (tanh — 

V a WHi \ a / CO sh^|^ 

a a 


a tanh 


rmeh 


p tanh^2^ 

a 


sin 


ritx 


J 


where r has even values. It Is evident that for all terms of this series In which r is 
greater than or equal to same number t, it 1 b approximately correct to write 


tanh 


cosh 


cosh 


rflPeh 

a 

rrtPeh 

a 

rrtcceh 

a 


tanh 


ritoeh 


a 

nt(Jeh 


-!• 


ritoeh 


-§• 


cs 


l 


i 


I 


BACA.nff Ho, i486 



Hence, approximately 


y a 



v sinh 
\ a ' coah 


ritmri 


-ft 


OMk ^ 

cl 


tanh 


r^hX 81 ”* ^ 


h \ 

y cosh 


rifBeb 


a tanh ^ - P tanh 


J 


Bln 


rnx 


+ x 


Q O 

a - ft) y . — 
r«t 


r«B 

■TT 


ntaeh / ritan rmiX 

V 


e a “ - e 


ritfteh / 

a [ a"“ a 

- e \© -e 


a ) 


sin 


ntx 


(46) 


where the summation is of even values of r and t is an even number and chosen so that 


tnfteh 


£ 2.4 


(49) 


A closer approximation results if t Is taken still larger than equation (49) wouli 
indicate. However, the approximation (equation (49)) appears to he reasonably good. 


After writing in exponential form the trigonometric factor sin occurring in 

the infinite series In equation (48), four infinite geometrical progressions are obtained, 
each of which is readily summed. It is found that 


H 

u> 
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T _ _Hi( 
j a 


rrtS 

a 


r= 

f 

e 


sinh einh Efa - 8 inh JS£ 8in h ^ 

— ] ain 223L 

, a Sinn £251^ 00ail £ 2 gih _ p Blnh rn£h coe}l rirgsh J a 


CslP t\ - e^ 1 aln(t - g)^ + e tp l [- a jja tX + e"^! Bin(t - 2)\] 
1 — 2e^l cos 2X + e^l 1 — 2e -2p l cos 2X + g*”^ 5 ! 


+ (a - p) 

[pin tX — e 2 ^ 2 alnft — 2)x] _ e tpg {j- gin tX + e~^ P ^ 8ln(t — 2)xl 


1 - 2s 272 cob 2X + o 7 '^ 


1 - 2e -2p2 cos 2X + e - ^ 32 


(50) 


H 

4=" 


vhere 


7i = f [a« (y - h) - 8] 

7 g => j [f* (r - h) - bJ 

p i * a |“ 6 (y + *0 + 8 ] 

p 2 “ ~ (y + h ) + 




(51) 


Usually, it is not necessary to take t larger than 12 or lU. However, in each oase 
equation (49) must he used in order to establish the value of t needed. 
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By applying this method to the other stress components and rewriting 
the expression for Xy, it results that equation (44) can "be written 

approximately as 


rrtS 


x .-?4 -| + ^ 2 


y - S 1 2 ' ^ D— ( P slDl1 rV 2 h C ° sh rV l 7 

— a sinh rv}h cosh rv gr ^ oos 

+ E^(3f(t,7i) + Pf(t, - Pi) -af(t,7 S ) - ctf(t, - r 2 )j L 


X x = - 


We 2 


rjrS 


T=c 


- ^. a ■ (a sinh 2V g h cosh tv^j 


rnx 


— P sinh rv<ih cosh rv^y ^ cos 
+ E £af(t,7i) + af(t, - Pi) - Pf(t, 7 2 ) - 0f(t, - P 2 )]J^ 


A 5 2 ) 


t— 2 


rit5 




^sinh rVgt- sinh rV-^y 


where 


\ rjtx 

- sinh rvih sinh rVgJ ) sln 
+ E[- if(t,7i) + if(t, - Pi) + if(t,7 2 ) “ ^ (4, - P2^j" 


J 


v -1 = 

1 a 


N, 


v 2 ' a 

* & ^ (t '" rwl) 

p(q.+i\)fT _ 2(q.-iX)l 




( 53 ) 
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i = Np-l; !)•]_, D 2 , y 1 , y 2 .^ P 3 ., P 2 , an d X are given by equations (45), (47), 

and (5l); and R means that the real part of the expression following It 
la to be taken. Correspondingly, equation (46) can be written approximately 
as 

STtS 


C 


f BJIO 

JP=£ e aT , 

A J> ^ — - ^3 cosh sv 2 h sinh sv^y 


1 


— a cosh sv^h sinh sVgy ^ cos 


1 


+ r{- (3f(m, n) + 0f(m, - P X ) + af(m,7 2 ) - af(m, - p 2 ) j 


** - "-H5 


fm—2 - 

l-'T — o 


Sjt8 


lira 


cosh sVgh sinh sv-j_y 


s?tx 

a 


i=T d 2 

— P cosh ev^h sinh sv^^ cos 
+ r[— of (rn^yj) + of(m, -P 3 _) + 0f(m,7 2 ) - 0f(m, - p 2 )J 

fxor-2 --•§2§. 

*y = ( COBh BV ^ COBh ^i 7 


r ( 54 ) 


— cosh sv-jh cosh 


svgy ) 


sin 


BTtX 


+ R j if(m, 7 X ) + if(m, - p^) - if (m, 7 g ) - if (m, - p 2 ) 

_ 0 4 J 

where s is odd; m «o when m is odd; Dp, Dg, y 72, Pp, P2, 

f, and i are given as 'before by equations (45), (47), (51), and C 53) J 
ami R means that the real part of the expression following it is to be 
taken. 


It is to be remembered that equation (52) for a small value of 5 
gives the stress components corresponding to the load distribution shown 
in figure 2(c), except that the loads are distributed over a 3mall area. 

The results for point loads are obtained by setting 8 equal to zero. 

(84 percent of the area under the curve of equation ( 31 ) is over a length 28 
on either side of x = 0.) Similarly, equation (54) corresponds to the load 
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■distribution shown in figure 2(b). Neither of these load distributions- 
alone would be of much practical importance. However, by adding the two 
distributions, the resulting stress distribution is that of an infinite 
beam loaded at equal intervals by equal concentrated forces acting in the 
upward end downward directions alternately, as shown by figure 2(a) . From 
the combined stress distribution, it is possible to draw some conclusions 
regarding beams of finite length subjected to concentrated loading. 


Reduction of Equations (52) and (54) for an Isotropic Beam 


The stress components given in equations (52) and (54) admit some 
simplification when they are applied to an isotropic beam. Although the 
isotropic case is not being considered in this report, it may be of 
interest to determine the forms to which equations ( 52 ) and (54) reduce in 
this case. For the isotropic case, P approaches a as a' approaches 1, 
and e also becomes equal to 1. For these values of a .and p, equa- 
tions (52) and (54) become indeterminate. By evaluating (noting that v, 

V 2 j 7} Vq) Pqj and P2 are functions of a and p), it results that for 
an isotropic beam, equation ( 52 ) becomes 



cosh ££ c osh =^-j cos a 


X 


+ *g (y - h)p(t, 7 ) - - (y + h)p(t,— p) - f (t, 7 ) - f(t, -p) 
fsinh (E3L sinh EEL + COB h &L) 



U 


V 


a J 


>(55) 


- CO sh E2* cosh EfL cos ES 
a a a j a 


+ E £ (y - h)P(t,7) - ^ (y + h)P(t,-p) + f(t,7) + f (t, -p) 


*y 


r?tS 

(ZSL alnh cosh ES2L 
8- * — 1 ' A An V a a 8- 


W 


|rb A 1 V 


— cosh sinh Ejp^sin 

Q, 8. J Q, 

* E [- X (y - h > p ( t ,5') - ^(rt 




J 
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where r Is even and 


t - ^ eTen ^ X = 7 = § “ h ) “ 8], p = * £(y + h) + &J 


p(q+iX) 

P(p ,4) = — [_P _ 2 (p - 2)e 2< l cos 


^1 — 2e 2< lcos 2X + 

+ (p — 4)e^ — (p + 2)e 2 ^~ 


2X 


+ 2pe ' * ' cos 2X — (p — 2)e 2 ^^ 


2(2q.-l\) 


f(p, q) = 


(a — 0) (l — 2e 2q - cos 2X + e 1 ^) 




56) 


- F 1 


A-, = + sirih cosh 

J_ cl cL a» 


J 


and B means that the real part of the expression following It Is to he 
taken. A closer approximation to the actual stress component In the 
Isotropic case results the larger the value of t Is taken. However, 

any value of t > appears from numerous computations to he reasonably 

ith 

satisfactory. Similarly, equation (55) becomes for an Isotropic beam: 
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7 a V;i a 2 


— slot 
a a 


^ sinh SISL | coe 
a a J a 





-p) 



4 


(57) 


where 


m > (odd) 

^ = _ B±nh e|h OOBh s|h 


s Is odd; 7, p, X, P, f, and i are given by equation (5^); and R means 
again that the real part of the expression following it is to "be used. 


EXPERIMENTAL VERIFICATION OF ANALYSIS 


T-n the vicinity of a concentrated load on a ‘beam, the distribution of 
the longitudinal shear over a cross section is markedly different from 
that predicted by the elementary . theory . An estimate of the accuracy of 
the me thod can be obtained by a comparison of an experimental determination 
of this distribution with that calculated by the mathematical method 
described herein. Such an experimental determination and a comparison 
were made and are described in the following sections. 
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Description of Test 

The plan was, broadly , to obtain a solid wood beam of uniform, struc- 
ture, to fix 4-7 metalectric gages on each £ace in the vicinity of the 
load point, to apply load and record strains by means of metalectric 
rosettes by using a 48— point recorder connected first to the gages on 
one face with one check gage on the other face, and then to repeat the 
application of load in Identical manner and record strains by using the 
recorder connected to the gages on the opposite faces. 

A clear, straight— grained Sitka spruce beam was selected. The rough 
piece was 3 inches by 10 inches by 16 feet- with the 10— Inch dimension in 
the tangential direction. The annual rings numbered about 18 to the Inch 
and their radius of— curvature was approximately 3 feet. The grain was 
almost parallel to the length of the beam throughout. One— half the rough 
piece was surfaced on four sides and trimmed on the ends to give a 
finished beam 1-99 by 9*37 by 96-03 inches. The specimen was then stored 
in a room of constant temperature and humidity for 3 weeks until its 
weight" became constant. 

It was desirable to obtain strains in three directions at definite 
points on the surfaces of the beam. The use of electric— resistance— type 
■strain gages appeared most feasible provided that- a short enough gage 
length can be obtained. Metalectric gages with g — Inch and ^4 — inch gage 

lengths and rosettes with 1— inch gage lengths are commercially available. 
Strains measured by these gages are average strains over the gage length, 
and, therefore, the 1— inch rosettes are not suitable for the purpose. 

Eosettes can be built up, with either the g — or inch gages, by mount- 
ing the gages on top of each other. The g— inch gages are l/4 inch wide, 

so that when they are superimposed on each othen, the first gage is 
shorter than the width of the gage beneath it - ; This was believed to be 

undesirable, and the g— inch gages were, therefore, not employed. 

The inch gages were l/8 inch wide and, therefore, this difficulty was 

not— encountered in their use. Rosettes were built up of these gages, and 
the effect of superimposing the individual gages on each other was found 
as follows. 

Three rosettes, built of three gages each, were mounted on the center 
line of a strip of clear Sitka spruce l/4 Inch thick, 1 inch wide, and 
about 24 inches long, which was then subjected to tension. The positions 
of the ‘gages are shown In figure 4, and the values of— the strains 
observed for a series of loads are given in table 1. In rosette' A, the 
longitudinal gage is applied directly on the wood. In rosette B, the 
longitudinal gage is between the other two gages. . In rosette C, the 
longitudinal gage Is superimposed on both of the other two gages. 
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Examination of the data, tabulated In table 1 shows that the strain 
readings were duplicated to within 0.00002 Inch per inch in successive 
loadings. The charts from the recorder were read accurately to about 
0.000008 inch per inch. Higher accuracy was not possible because of the 
thickness of the recording trace. When the traces overlapped prior 
traces, the readings were less accurate. 

The data also indicate that the position of the gage in the roBette 
does not significantly influence the strain recorded. The longitudinal 
strain measured by gage 1 in rosette A (fig. 4) can be assumed to be 
correct because the gage was mounted directly on the wood. Strains 
measured by gage 2 of rosette B, which was superimposed on one other 
gage, showed increases , compared with gage 1 of 0.000010 inch per inch in 
the first and 0.000016 inch per inch in the second loading at a load of 
750 pounds. The longitudinal strains measured by gage 3 of rosette C at 
the same successive loads were 0.000060 and 0.000068 inch per inch 
greater than those measured by gage 1. Thus these data indicate that the 
superimposed longitudinal gages undergo greater strains than the gage 
mounted directly on the wood. However, the gages measuring strain at 4^ 
to the longitudinal direction do not confirm thiB indication. Gage 6 was 
applied directly on the wood. Gage 4, which was applied on one other 
gage, yielded values of strain 0.000040 and 0.000050 inch per inch less 
than those of gage 6 ; and gage 5 j which was superimposed on two other 
gages, yielded values 0.000050 and 0.000040 inch per inch lesB than those 
of gage 6. Gage 8 was also applied directly on the wood. Gage 9, which 
was superimposed on one other gage, yielded values 0.000060 inch per inch 
less than that of gage 8; and gage f, which was superimposed on two other 
gages, yielded values 0.000030 inch per inch greater than the value of 
gage 8. These results indicate that the position of a gage (bottom, 
center, or top) in a built-up rosette apparently has no consistent effect 
on the recorded strains. Furthermore, the differences in the magnitudes 
of the strains read by gages oriented in the same direction are so small 
that they might very well be differences in actual strain from point to 
point in the specimen. It is assumed, therefore, that rosettes of this 
type yield values of strain that are sufficiently accurate for the purpose 
of this report. 

Gages of this type were mounted on the beam while it was approaching 
its equilibrium, moisture content in the humidity room. Their positions 
are shown by the sketch (fig. 5) and by the coordinates in table 2. The 
rosettes were built up with the gage measuring strain in the longitudinal 
direction of the beam applied directly to the wood and centered on the 
point, shown in figure 5* The gage measuring strain at 90° to the longi- 
tudinal direction was superimposed on the first., gage and centered on the 
same point. The third gage, which measured strain at 45° to the longi- 
tudinal, was superimposed on the first two gages and centered on the 
same point.' Additional gages with 1— inch gage length were applied to 
measure . longitudinal strains in the vicinity of the. neutral axis as 
shown. ‘‘ 
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Point-s symmetrically arranged at out the line of action of the load 
were accurately laid out— on both faces of the beam by using a square 
and a scale graduated to 0.01 inch. Application of the gages required 
that a prime coating of glue be allowed to dry on the wood. Fresh glue 
was spread on both the gage and the beam, and then the gage was firmly __ 
pressed into position, the squeezed— out glue was removed, and a weighted 
piece of sponge rubber to apply pressure was left on the gage for 48 hours. 
The layout lines were covered by the first gage applied, so that the lines 
were redrawn on top of the gage for applying the superimposed gage. How- 
ever, after all gages were applied and the coordinates were measured for 
the center of each gage, it was observed that'some individual gages were 
centered as much as 0.04 inch from the point desired. In comparison with 
the dimensions of the beam, the location of gages was considered satis- 
factory and was within 0.03 inch of the average values given in table 2. 

The beam was tested in the room of controlled temperature and 
humidity after equilibrium moisture content was attained. A four— screw 
mechanical testing machine of 10, 000-pound capacity waB used with a static- 
bending Jig centered on the weighing platform. The jig was made of a 
pair of ip— inch I— beams bolted together with spacers to give a £ — inch 

clearance between flanges. Laterally adjustable supports were spaced 
symmetrically from the center of-^bhe jig to give a 6-foot span and were_._ 
bolted to the top flanges. The beam was then centered on the supports, 
and roller-bearing plates were inserted between the beam and the laterally 
adjustable supports . Load was applied to the center of the beam by means 
of a spherical head and a hard-maple bearing block cut in the shape of a 
cylindrical segment. 

A great amount of difficulty was encountered in attempts to apply 
the load. Several loading blocks were each tried in several preliminary 
loadings, butrthe observed values of strain were far from symmetrical 
about the plane in which load was applied. The slightest~change in the 
position of the loading block produced large change In the strains on the 
faces of the beam. Consecutive loadings under conditions reproduced as 
closely as possible for each load did not "produce identical strains. The 
plan to observe strains in each face of the beam in two consecutive load- 
ings was therefore discarded. The 1— Inch gages in the vicinity of the 
neutral axis could not be read with sufficient accuracy to yield useful 
data; 48 gages were, therefore, selected for measurement of strain during 
a single application of the load. These gages are shown in circles in 
figure 5- 

Before the final data were taken, a jointer cut 0.05 inch deep was 
taken along the top of- the beam to remove any material that could have 
been overstressed beneath the load block. during the previous .applications 
of load. A new hard-maple load block was 'turned. In a wood lathe to an 
11— inch radius to insure a true cylindrical segment. The beam was again 
set up and small loads, less than 200 pounds, were applied to check. the 
centering of the loading block against the sy mm etry of the gage readings. 
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Movements of 0.01 Inch of the loading "block were sufficient to produce 
definitely nonsymme.trlcal distributions of strain. When the -loading 
block. was centered to . give approximately symmetrical strains, - the beam 
was tested as follows. 

The 48— point recorder was' adjusted for the initial reading of each 
gage while a load of 15 pounds was maintained on the beam.' " Strains were 
then recorded for gages on both faces of the beam while consecutive loads 
of 100, 200, 300, 400, 500 , and 600 pounds were maintained. 

Observed strains were recorded while a constant deflection was main- 
tained on the beam. Because the load was comparatively small and Only "a 
few minutes were required to record the strains, plastic flow "of the 
material of the beam probably did not influence the strains observed. 
Certain gages registered such small increments of strain (less than 
0.000008 in. /in.) that the recorded traces overlapped and made accurate 
interpolation of the strain difficult . The readings taken for small 
increments of strain are, therefore, not so reliable as those taken for 
large increments. 

After the beam was tested, specimens for the determination of its 
elastic properties were cut from it. They were taken from the central 
part of the beam. For the plate— shear test, these specimens consisted of 
two specimens measuring l/4 by 9^- by inches each; and for the com- 
pression test, the specimens consisted of two specimens parallel to the 
grain and measuring 2 by 2 by 8 inches each and two specimens perpendicular 
to the grain and parallel to the depth of the beam and measuring 2 by 2 
by 8 inches each. 

From the static bending tests the modulus of elasticity only was 
determined. The modulus of rigidity in planes parallel to the faces of 
the beams was dete rmin ed by the plate-shear tests. Modulus of elastic- 
ity in compression and Poisson’ s ratios in the TE— , TL-, LT— , and re- 
directions were obtained by the compression tests. 


Presentation of Data - - 

Data obtained from the test of the beam are recorded in table 3* In- 
column (l) are listed the west— and east— face gages shown in figure 5* The 
gage readings at 15 , 100, 200, 300, 400, 500, and 600 pounds are given in 
columns ( 2 ) to (8), respectively. Dimensions of the beam at time of test 
were 1.99 by 9*32 by 96 . 03 inches. The strains tabulated were read and 
checked from 1 the chart plotted by the 48— point strain recorder. 

.Readings of each gage are plotted in figure 6 and are grouped in 
accordance with distance from the plane in which the load was applied and 
with distance from the center line of the beam. The slope of these curves' 
from 0 to 500 pounds was used to determine the strains" tabulated ' in 
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column ( 9 ) of table 3 • From these strains the shear strains in the sy— plane 
were computed for the points at which rosettes were located. In table 2 
the x— and- y— coordinates of the gages are given in columns 2 , 3 , 6 , and 7 , 
and the shear strains are tabulated in columns (4) and ( 8 ) . The strains In 
both faces and at points 0.4 inch on each side of- the load were averaged 
for each group of gage positions at symmetrical locations. These average 
values are given at the end of table 2 . ' ~ 


Shear strains in the xy— plane of the beam when subjected to a 500— pound 
load were computed by the mathematical method previously presented for 
points 0.4 Inch from the plane in which load was applied. Computations 
were based on the elastic properties given in table 4 obtained by tests of 
coupons cut from the beam after test. The strains were calc ula ted for 
three distributions of load (6 equal to 0 , 0 . 125 , and 0.250 inch) by 
obtaining the shear stresses 2 y from the third par 1 r~of equations ( 52 ) 
and (54) and dividing their sum by the modulus of rigidity. This method 
of obtaining shear strains from shear stresses Is valid only when the 
strains associated with the grain direction and the radial or tangential 
direction are required. (See references 8 and 9*) 


For example, the strain computed at x = 0.40 and y = 1.48 with 8=0 
was obtained as follows. First the value of_ |- was found from equation (4) 
to be 0-03594, and from equation (49) the value of t was established 
as 8 . Then from equation (52) the term 

2 rrt& 

. Q - y 8 ! - ^sinh rv^h sinh rV^y — sinh rv^ sinh rvgy^ sin 

was evaluated and found to be 0.001002. The real part of the term 
- — E^-if (t, 7l ) + if (t, -Pi) + if(t, 7 2 ) - lf(t, -P 2 fj 


of equation (52) was found to be 0.002806. By adding these terms and 
multiplying by W the stress X^ was found equal to 0.003808W. In a 
similar manner for equation (54), m was first-nerstablished equal to 7 . 
The term 


Srt5 


fr / — — - (cosh sVph cosh ev-,7 — cosh SV h cosh sv 0 y'\ sin 

s=l d 2 \ x 2 / 


srtx 


and the term 


|-Ejif(m, 7q) + if(m, -p x ) - if(m,7 2 ) - if(m, -p 2 )J 


were found to be 0.004437 and .0.003925, respectively. On summing and 
multiplying by W, the value of— X y from equation (54) was found equal 
to O.OO 8362 W. The sum of the stresses computed by equations ( 52 ) and (54) 
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gave the total shear stress X_ at the point equal to 0.01217W. By 
substituting the lead per inch? of width (500 lb/ 1.99 in.) for W and 
dividing by the modulus of rigidity * 1 ™., the strain e r , was found 
equal to 33-72 X 10 -6 . ^ 

Computed strains for several points are given in table 5 in which 
columns (l) and ( 2 ) give the x— and y— coordinates of the points and 
columns (3)* (4), and (5) give the strains, respectively, with S equal 
to 0 , 0 . 125 , and 0.250 inch. 

A comparison of the computed and the observed distribution of shear 
strains is presented in figure 7 in which three curves, one for each 
value of 8 , and the average observed strains are plotted with distance 
from longitudinal center line as ordinate and the strains at points 
0.4 inch from the plane of loading as abscissa. 


RESULTS 


The results of the computations given in table 5 and the average 
results of the tests given in table 2 are plotted in figure 7- In this 
figure the ordinates are distances upward from the center line of the 
beam and the abscissas are shear strain. The distribution of strain 
obtained by the elementary theory is also shown. The difference between 
the two methods is marked; that is, the method developed herein exhibits 
a high— stress concentration near the top of the beam and the elementary 
method yields a low maximum at the center of the beam. 

The data from the test exhibit a stress concentration similar to 
that obtained by the more accurate theory. The measured strains at the 
two points 1.48 and 2.72 inches from the center of the beam agree with 
the theory within the accuracy of the experiment. The strains at these 
two points are so small that the corresponding traces on the chart of the 
recorder overlapped and could not be accurately read. The strain at 
3.72 inches from the center of the beam agrees very well with the curve 
for 8=0 that is for a concentrated load. Of course, the load was 
not concentrated in the test but was applied over a length of about 
0.4 inch by the cylindrical load block. However, at a distance from the 
region over which the load is applied the difference in the effect of a 
truly concentrated load and one applied over a small region should be small. 
The strain at 4.48 inches from the center of the beam shows the effect cf 
the distribution of the load. At this point, the theory approximately 
agrees with experiments if 8 is given the value of 0.250 inch. The 
strain is considerably less than that due to a concentrated load. 

If the actual distribution of the load on the beam were known, a more 
accurate solution for the strain distribution near the load could be 
obtained by integrating the solution for a concentrated load of varying 
intensity over the loaded part of the beam. However, this distribution 
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is difficult to determine and usually is not- known, and, therefore, further 
refinement of the method seems futile. 


CONCLUDING- REMARK 


It appears that the shear strains computed hy a derived mathematical 
method were verified, as closely as could "be expected, by experiment~on 
a wood beam of-rectangular cross section. 


Forest Products Laboratory, Forest Service 
U. S. Department of Agriculture 

Madison, Wis., February 18, 194-7 
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TABLE 1.- STRAIKS OBSERVED DT EOSETTBS MOUHTED OH A STOOP OF WOOD 


SUBJECTED TO A TEH3XLK LOAD 


Load. 

(ll) 




Strain 
(la. /in.) 

- 




Gage 1 

■ 

Gage 2 

— 
Gage 3 

Gage 1 

Gage 5 

Gage 6 

Gage 7 

Gage 8 

Gage 9 

Seoond ^plioaticm of load 

50 

100 

150 

200 

250 

300 

350 

1*00 

1+50 

500 

550 

600 

650 

700 

750 

108 x ict 6 
216 
328 
10(8 
561 
680 
800 
901 

1032 

1118 

1268 

1392 

1516 

16I0 

0 , 

100 x 10“° 
212 
328 
10*8 
561 
681 
796 

■ 916 
1036 
1160 
1280 
lloo 
1521 
1656 

0 

112 X 10 "° 
22l 
336 
156 
576 
700 
012 
936 
1060 
1181 
1312 

lllo 

1518 

1708 

° -6 
20 X 10 "° 

1(0 

50 

60 

80 

100 

120 

130 

160 

170 

190 

200 

220 

210 

10 x 10 -6 
30 
50 
60 
80 
100 
120 
130 
160 
180 
200 
210 
230 
250 

0 * 
10 x 10“° 

30 

50 

60 

90 

110 

130 

llO 

170 

190 

220 

2lO 

270 

290 

° 

60 X 10 “° 
120 
170 
210 
270 
320 
370 
110 
l70 
510 
560 
600 
650 
TOO 

60 x ior6 
110 
160 
200 
250 
300 
3lO 
390 
I3O 
1(80 
530 
570 
530 
670 

50 x lo" 6 
100 
150 
180 
230 
270 
300 
350 
390 

I30 

180 

520 

57O 

610 

First application of load 

50 

0 

0 

0 

0 

0 

0 

0 

0 

0 

100 

— 

100 

100 

20 

10 

10 

lo 

lo 

10 

150 

— - 

210 

220 

10 

20 

20 

100 

90 

90 

200 


320 

330 

50 

10 

50 

150 

150 

llo 


— 

150 

150 

80 

60 

60 

210 

190 

180 


— 

570 

570 

80 

80 

90 

260 

250 

230 

350 

— 

690 

690 

100 

100 

no 

300 

300 

270 

1(00 


810 

810 

120 

120 

130 

350 

3 ^ 

310 

1*50 

— 

910 

93O 

150 

130 

150 

120 

390 

350 

500 

1030 

1030 

1060 

170 

150 

170 

170 

I30 

390 

550 

II60 

1170 

1180 

190 

170 

190 

520 

I70 

130 

600 

1270 

1290 

1310 

200 

190 

210 

560 

520 

180 

650 

iloo 

llio 

lHo 

220 

210 

2lO 

610 

560 

520 

700 

1520 

1520 

1560 


230 

260 

650 

610 

560 

750 

i 64 o 

1650 

1700 

mmi 

250 

300 

TOO 

670 

1 

610 
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TABES 2,- OBSERVED STRAINS, AVERAGE VALUES, AND COORDINATES EOR POINTS AT 
WHICH GAGES WERE LOCATED ON BEAM TESTED (AXES SHOWN IN FIG. l) 


Gage 

Coordinates ' 

Shear 

strain, 

© 

( in. /in.) 

Gage 

Coordinates 

Shear 

strain, 

e sy 

(in. /in.) 

(in.) 

7 

(in.) 

(to.) 

y 

(in.) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


West 

face 



East 

Face 


1, 2, 3 

0.97 

4.48 

442 x 10" 6 - 

1, ?, 3 

-1.03 

4.48 

-426 x 10-6 . 

4, 5, 6 

.38 

4.49 

1140 

4, 5, 6 

-.40 

4.43 

— l4oo 

7, 8, 9 

-.40 

4.48 

-914 

7, 8, 9 

.40 

4.47 

946 

10, 11, 12 

-l.oi 

4.48 

-314 

10, 11, 12 

-.40 

3.72 

-316 

13, 14, 15 

• 37 

3.71 

376 

13, 14, 15 

.41 

3.73 

■338 

16, 17, 18 

-.40 

. 3.72 

-260 

16, 17, 18 

.40 

2.72 

— — 

19, 20, 21 

• 37 

2.72 

50 

19, 20, 21 

-.40 

1.48 

-50 

"22, 23 , 24 

-.40 

2.73 

-58 

22 

.99 

• 23 

— 

25, 26, 27 

.38 

1.48 

— 

23 

.99 

.00 

— 

28, 29, 30 

-.40 

1.48 

-36 

24 

1.00 

-.19 

ma am mm 

31 

5.01 

.02 

— 

25 

.01 

.22 

— — 

32 

5.01 

.21 

---- 

26 

-.01 

-.17 

— 

33 

5.00 

-.17 


27 

.01 

.39 

---- 

34 

.94 

.22 

— 

28 

.00 

.01 

— — 

35 

• 95 

-.17 


29 

-1.04 

.21 

— 

36 

.94 

.02 

— 

30 

-1.03 

.00 


37 

—•^3 

.23- 

— 

31 

-1.04 

-.17 

— 

38 

-.03 

-.17 

— 

32_ 

.00 

4.48 


39 

-.04 

.38 

— 

Top of ‘beam 

| 1 ^| 

— 

40 

-.04 

.03 

— 





4l 

-1.04 

-.17 






42 

-1.05 

.23 

— 





43 

■ -1.05 

.03 

— 

Average 

±0.40 

4.48 

1100 

44 

-5.00 

.23 

— — 

values 

± .40 

3.72 

322 

45 

-4.99 

.02 

— 

("both 

± .40 

2.72 

54 

46 

-5.00 

-.17 

— 

faces) 

± .40 

1.48 

43 

47 

-.03 

4.48 

— 
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TABLE 3.— LQAL-SERAIN DATA OBTAINED IN URAL TEST OT BEAM 


Oaga 

Gage readings, 
(in. /in.) 

Strain Incre- 
ment from 
0 to 500 lb 
(in. /in.) 

15-lb load 

100-lb load 

200-lb load 

30O-lb load 

400-lb load 

500-lb load 

600-lb load 

JiL 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

1ST 

582 x 20-6 

560 X 20 -6 

536 X 20^ 

505 x 20 -6 

479 x 20 r-6 

452 x 20-8 

432 X 20 -6 

-272 x 10-6 

2V 

339 

340 

341 

341 

342 

342 

342 

6 

3 w 

339 

.341 

350 

360 

370 

381 

384 

88 

4tf 

585 

570 

553 

537 

520 

505 

499 

-170 

5W 

935 

896 

832 

767 

700 

625 

553 

-670 

6v 

632 

642 

655 

673 

686 

700 

699 

150 

7V 

550 

544 

534 

528 

527 

528 

532 

-e4 

0W 

960 

934 

878 

810 

720 

615 

495 

-748 

9W 

646 

660 

675 

685 

690 

698 

698 

76 

1CW 

780 

760 

725 

702 

678 

648 

615 

-274 

11V 

575 

577 

577 

575 

575 

573 

573 

0 

I2V 

427 

428 

429 

'430 

430 

438 

441 

20 

13V 

587 

580 

573 

570 

566 

562 

554 

-40 

14V 

1000 

970 

915 

863 

812 

763 

720 

—696 

15V 

>47 

735 

718 

699 

681 

661 

647 

-180 

l6v 

682 

678 

671 

667 

660 

654 

648 

-52 

17V 

792 

761 

707 

660 

610 

540 

520 

-480 

l8w 

735 

728 

710 

700 

686 

672 

660 

-136 

1SW 

684 

679 

672 

664 

660 

652 

647 

-64 

sow 

700 

685 

66 2 

640 

620 

600 

580 

-210 

21V 

642 

635 

620 

612 

600 

589 

578 

-112 

22W 

582 

578 

568 

562 

556 

550 

542 

-70 

23W 

908 

892 

862 

840 

818 

790 

770 

-244 

2 4W 

840 

836 

815 

805 

790 

778. 

762 

-128 

-19E 

590 

588 

585 

583 

582 

581 

581 

-18 

20® 

635 

616 

598 

580 

565 

578 

532 

-160 

21E 

570 

565 

558 

553 

547 

553 

532 

-64 

28W 

692 

690 

685 

680 

678 

672 

670 

-44 

2SW 

512 

502 

490 

480 

464 

452 

440 

-128 

30V 

654 

650 

642 

636 

630 

622 

616 

-68 

11E 

755 

6t3 

625 

570 

520 

465 

4l8 

-538 

12E 

726 

708 

700 

688 

678 

665 

658 

-120 

lfiE 

695 

693 

690 

687 

683 

680 

679 

-60 

7® 

582 

554 

525 

500 

472 

448 

423 

-254 

8E 

790 

748 

720 

697 

668 

635 

595 

-272 

9E 

452 

487 

511 

531 

557 

578 

602 

210 

4E 

578 

545 

522 

500 

480 

458 

442 

-220 

6e 

898 

842 

803 

760 

710 

660 

600 

-460 

5E 

270 

312 

342 

383 

421 

455 

481 

360 

IE 

585 

559 

533 

500 

478 

450 

425 

-270 

ZE 

405 

4o4 

402 

402 

401 

400 

399 

-16 

3® 

351 

360 

365 

372 

380 

385 

392 

70 

13E 

591 

592 

593 

593 

593 

593 

593 

4 

ite 

890 

807 

745 

687 

630 

573 

520 

-574 

15® 

540 

526 

524 

502 

493 

482 

472 

-116 

10 E 

695 

694 

690 

688 

686 

682 

680 

-28 

32® 

692 

720 

588 

455 

327 

211 

111 


1+7W 

1003 

875 

682 

505 

337 

210 

100 







TABLE 4.- ELASTIC PROPERTIES OF SITKA SPRUCE BEAM 


FROM TESTS OF COUPOHS 


*rv .1 

Test value 

J 

Average 

value 

Property - 1 

Coupon 1 

Coupon 2 

Modulus of rigidity In shear (lb 
per sq In.) 

90,560 

90,820 

90,690 

Modulus of elasticity In compres- 
sion 

E^ (lb per sq In.) 

i, 6 ia, 000 

1 , 757,000 

1 , 699,000 

Erp (lb per sq in.) 

34,870 

40,950 

37,910 

Poisson's ratio 
a LT 

0.b6 

0.50 

0.L8 

°IR 

.21 

.25 

.23 

0 TL 

.01 

.02 

.02 

°TR 

.66 

.66 

.66 


^L, T, and R refer to direction of grain, direction tangent to growth 
rings, and direction perpendicular to growth rings, respectively. 


oo 
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TABIE 5.- COMPUTED STRAINS IN BEAM AT 500-P0UND LOAD FOR 
6=0, 0.125, AND 0.25 AND COORDINATES OF POINTS 


Coordinates 



Computed strain, e 
( in. /in. ) 

xy 

6=0 In. 

6 = 0:125 in. 

5 = 0.250 In. 

33-72 x 10-6 

33-25 x 10“6 

30.87 x 10“6 

80.46 

74.39 

68.95 

146.22 

131.15 

118.17 

312.10 

266.53 

229.38 

777.67 

61O.7O 

488.07 

1490.50 

1088.34 

813.02 

2591-00 

1832.24 

1280.66 


1922.92 

1344.37 


1825.47 

1308.34 


1361.33 

1059.90 



- .06 

.06 
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Figure 4.- Gage positions in built-up rosettes 
on tension specimens for check test. 
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Figure 5.- Positions of electric -resistance -type strain gages on. 
both faces of beam. (See table 2 for x- and y-coordinates of 
gage centers.) 





y-coordinate, in. 



Figure 7.“ Comparison of observed strains with curves showing strains computed 

- tm mm i /\ n . A 1. _ A GO mJ. CAH nAn+AH IaO TJrvtti 

al on g me line x = u.^u irom y = u w y = ai uvwtwuiw ^mcj. iucw, 

the elementary theory and the advanced mathematical theory of this report were 

used, with 6 equal to 0.0, 0.125, and 0.250 inch. 
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